Higher-order effects and ultra-short solitons in left-handed metamaterials 
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Starting from Maxwell's equations, we use the reductive perturbation method to derive a second- 
order and a third-order nonlinear Schrodinger equation, describing ultra-short solitons in nonlinear 
left-handed metamaterials. We find necessary conditions and derive exact bright and dark soliton 
solutions of these equations for the electric and magnetic field envelopes. 



Electromagnetic (EM) properties of metamaterials 
with simultaneously negative permittivity e and perme- 
ability \i have recently become a subject of intense re- 
search activity. Such metamaterials were experimentally 
realized recently in the microwave regime, by means of 
periodic arrays of small metallic wires and split-ring res- 
onators (SRRs) [l| . Many aspects of this class and other 
related types of metamaterials have been investigated, 
and various potential applications have been proposed 
0. So far, metamaterials have been mainly studied in 
the linear regime, where e and /i do not depend on the EM 
field intensities. Nevertheless, nonlinear metamaterials, 
which may be created by embeddingan array of wires 
and SRRs into a nonlinear dielectric [1, 0, Q , may prove 
useful in various applications. These include "switching" 
the material properties from left- to right-handed and 
back, tunable structures with intensity-controlled trans- 
mission, negative refraction photonic crystals, and so on. 

EM wave propagation in nonlinear metamaterials can 
be described by two coupled nonlinear Schrodinger (NLS) 
equations for the EM field envelopes [6|. Thus, bright- 
bright and dark-dark vector solitons of the Manakov 
type [7|] are supported in the right-handed (RH) and 
left-handed (LH) regimes, respectively @. These find- 
ings paved the way for relevant studies, e.g., modula- 
tional instability Q, and bright-dark vector solitons @ 
in negative-index media. Additionally, a scalar higher- 
order NLS (HNLS) equation was derived in [10( (as- 
suming nonlinear response only in the electric proper- 
ties of the metamaterial), and was subsequently studied 
[TTI . \v\ . Coupled HNLS equations were also derived E3 , 
where higher-order dispersion and nonlinear effects were 
included. However, the relative importance of these ef- 
fects was not studied in Ref . EH , although such an inves- 
tigation should provide the necessary conditions for the 
formation of few-cycle pulses in nonlinear metamaterials. 

In this work, we present a systematic derivation of NLS 
and HNLS equations for the EM field envelopes, as well as 
ultra-short solitons for left-handed (LH) metamaterials. 
In p articular, we use the reductive perturbation method 
[14| to derive from Faraday's and Ampere's Laws a hier- 
archy of equations. Using such an approach, i.e., directly 
analyzing Maxwell's equations, we show that the electric 



field envelope is proportional to the magnetic field one 
(their ratio being the linear wave- impedance) . Thus, for 
each of the EM wave components we derive a single NLS 
(for moderate pulse widths) or a single HNLS equation 
(for ultra-short pulse widths), rather than a system of 
coupled NLS equations (as in Refs. [E [H, [E [HI). The 
HNLS equation, which incorporates higher-order disper- 
sive and nonlinear terms, generalizes the one describ- 
ing short pulse propagation in nonlinear optical fibers 
EE EE El El- Analyzing the NLS and HNLS equa- 
tions, we find necessary conditions for the formation of 
bright or dark solitons in the LH regime, and derive an- 
alytically approximate ultra-short solitons in nonlinear 
metamaterials. 

We consider lossless nonlinear metamaterials, charac- 
terized by the effective permittivity and permeability 0] , 
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where to p is the plasma frequency, F is the filling factor, 
^onl is the nonlinear resonant SRR frequency [3j, while 
E and H are the electric and magnetic field intensities, 
respectively. In the linear limit, er> — > 1 and luqnl - > 
L0 Tes (where w rcs is the linear resonant SRR frequency), 
and LH behavior occurs in the frequency band Lo rcs < 
to < min{(xv, 



p ,lo m }, with to M = w ros /Vl - F, provided 
On the other hand, a weakly nonlinear 
behavior of the metamaterial can be approximated by 
the decompositions [E H, EE EH : 



that LOp > to- 



e(to) = e L {Lo) + e NL (uj;\F,\ 2 ), 
fj,(co) = fi L (uj) + fi NL {uj;\tl\ 2 ) 



(3) 
(4) 

»)]. 



where e L = e (l - uj 2 /lo 2 ), [i L = fi [l - Flo 2 / (lo 2 - 
while the nonlinear parts of the permittivity and perme- 
ability are given by [EIEEE|: ^l(|E| 2 ) = e a|E| 2 , 
and ^jvl(|H| 2 ) = /i /3|H| 2 ; here, a = ±E~ 2 and $ are 
the Kerr coefficients for the electric and magnetic fields, 
respectively, E c being a characteristic electric field value. 
The approximations ©-{H are physically justified con- 
sidering that the slits of the SRRs are filled with a non- 
linear dielectric (EH. Generally, both cases of focusing 
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FIG. 1: (Color online) The linear parts of the relative mag- 
netic permeability, hl/ho [solid (red) line], and the electric 
permittivity, el / eo [dashed (blue) line] as functions of fre- 
quency, for F — 0.4 and lo p — 2ty x 10 GHz. 

and defocusing dielectrics (corresponding, respectively, 
to a > and a < 0) are possible. The magnetic Kerr 
coefficient (3 can be found via the dependence of [i on the 
magnetic field intensity 0, Q. Here, fixing F = 0.4 and 
Lo p = 2ir x 10 GHz, we will perform our analysis in the 
frequency band from 2tt x 1.45 GHz to 2tt x 1.87 GHz, 
considered also in Ref. L 6'j. In this band, SRRs are LH 
media (with cl < and ttj, < - see Fig. 1), a may be 
either positive or negative, while (3 is positive @. Notice 
that, in principle, cnl and ^jlnl may depend on both in- 
tensities |E| 2 and |H| 2 ; such a case can also be studied 
via the analytical approach we use below. 

We consider the propagation along the +z direction 
of a x- (y-) polarized electric (magnetic) field, namely, 
E(z,t) = xE(z,t) and H(z,t) = Y H (z,t). Then, using 
the constitutive relations (in frequency domain) D=eE 
and B=/iH (D and B are the electric flux density and 
the magnetic induction), Faraday's and Ampere's Laws 
respectively read (in the time domain): 



d z E=-d t ((i*H), d z H = -dt{e*E), 



(5) 



where * denotes the convolution integral, i.e., f(t)*g(t) — 

J-°° f( T )g{t ~ t)cLt. Note that Eqs. (JSJ) may be used in 
either the RH or the LH regime: once the dispersion rela- 
tion ko = fco(^o) (for the wavenumber fcrj and frequency 
u>o) and the evolution equations for the fields E and H are 
found, then ko > (fcrj < 0) corresponds to the RH (LH) 
regime. Alternatively, for fixed fco > 0, one should shift 
the fields as [E, H] T — > [±E, TH] T , thus inverting the 
orientation of the magnetic field and associated Poynt- 
ing vector. Here, we will assume that the wavenumber 
ko [see Eq. ([16]) below] will be ko < for the LH regime. 

Now, we consider that the fields are expressed as 
[E(z,t),H{z,t)] T = [q{z,t),p(z 7 t)} T exp[i(koz - w t)], 
where q and p are unknown field envelopes. Nonlinear 
evolution equations for the latter can be found by the 
reductive perturbation method [l4[ as follows. First, we 
assume that the temporal spectral width of the nonlinear 
term with respect to that of the quasi-plane- wave disper- 
sion relation is characterized by the small parameter e 
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Then, we introduce the slow variables: 
Z = e 2 z, T = e(t-k' z), (6) 



where k' = v~ is the inverse of the group velocity (here- 
after, primes will denote derivatives with respect to too)- 
Additionally, we express q and p as asymptotic expan- 
sions in terms of the parameter e, 

q{Z,T)=qo{Z,T)+sq x {Z,T)+s 2 q 2 {Z,T) + --- , (7) 
p(Z,T) = Po (Z,T) + ep! {Z,T) + e 2 p 2 {Z,T)-\ , (8) 

and assume that the Kerr coefficients a and (3 are of order 
0(e 2 ) (see, e.g., (3, EE El)- Substituting Eqs. ©-© 
into Eqs. ©, using Eqs. Q, ((4]), and ([6]), and Taylor 
expanding the functions ej,, and p,L, we arrive at the 
following equations at various orders of e: 
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C(e 3 ) : Wx 3 = -iW'd T X2 + -W"fi£xi 
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Axi + iBxq, 
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with Xi= [qi,pi] (i = 0, 1, 2,3) unknown vectors, and 
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with -k denoting complex conjugate. To proceed further, 
we note that the compatibility conditions required for 
Eqs. ([5|)-(fl2l to be solvable, known also as Fredholm 
alternatives [l4l [l6j . are LWx^ = 0, where L = [1,Zl] 
is a left eigenv ector of of W, such that LW = 0, with 
Zjj = \J P-l/cl being the linear wave-impedance. 

The leading-order Eq. ((9]) provides the following re- 
sults. First, the solution xq of Eq. J9j has the form: 



x =R.<f>(Z,T), 



(15) 



0. 



where <f>(Z, T) is an unknown scalar field and R 
[1, Z^ 1 ] T is a right eigenvector of W, such that WR = 
Second, by using the compatibility condition LWx = 
and Eq. (fl5)l . we obtain the equation LWR = 0, which 
is actually the linear dispersion relation, 
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(ez, and m, are evaluated at cjq). Note that Eq. (fl6|) 
is also obtained by imposing the nontrivial solution con- 
dition detW = 0. Third, the EM field envelopes are 
proportional to each other, i.e., qo = Po^l- 



3 




10.6 11 
00 = 2ti f (GHz) 



TABLE I: Conditions for the formation of bright or dark soli- 
tons (BS or DS) for the NLS Eq. QTQ. 
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FIG. 2: (Color online) The GVD coefficient k" = d 2 k/du 2 
as a function of frequency in the left-handed regime. 

At C^e 1 ), the compatibility condition for Eq. (|10[) re- 
sults in LW'R = 0, written equivalently as: 



2k k' = u 2 (e L f/ L + e' L ^ L ) + 2uj e L fj, L . 



(17) 



This is actually the definition of the group velocity 
V g = 1/A;q, as can also be found by differentiating Eq. 
(fT6|) with respect to u. Furthermore, using Eq. (fT5|l . Eq. 
03) suggests that the unknown vector X! has the form, 



xi = iR'd T <l>(Z, T) + Rip(Z, T), 



(18) 



where ip(Z,T) is an unknown scalar field. 

Next, at order 0(e 2 ), the compatibility condition for 
Eq. QI]), combined with Eqs. 03 and ([15]). yields the 
following NLS equation, 



0, 



(19) 



where kg is the group-velocity dispersion (GVD) coef- 
ficient, as can be evaluated by differentiating k' in Eq. 
OH) , and 7 = (wg/2fc )(e aML + fi (3e L Z^ 2 ). Note that 
once </> is obtained from the NLS Eq. ([T9]). the EM field 
envelopes are determined as <?o = 4> an d Po = Z^cf) [see 
Eq. (|15p ]. similarly to the case of a linear medium. 

Finally, to order 0(e 3 ), we use the compatibility con- 
dition for Eq. 02) i as weu as E q s - 03), 03) and 03, 
and obtain a NLS equation, incorporating higher-order 
dispersive and nonlinear terms. This equation describes 
the evolution of ip, and yet contains <f>, which in turn 
obeys Eq. (fl~9|) . Instead of considering this system of 
two equations, we follow [HI, EH, HH and introduce a new 
combined function <f> = (j) + eip. This way, combining the 
NLS equations obtained at orders 0(e 2 ) and 0(e 3 ), we 
find that $ obeys the HNLS equation: 



is 
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(20) 



For e = 0, the HNLS Eq. flU is reduced to the NLS 
Eq. 02), while for e ^ generalizes the higher-order 



NLS equation descri bing ultra-short pulse propagation 
in optical fibers [TH, [161 [TtT. Il8j (where dispersion and 
nonlinearity appear solely in the dielectric properties). 
As in the NLS Eq. 03), E q. © provides the field $ 
which, in turn, determines the EM fields at order 0(e 3 ) 
as qo+eqi = <f> andpo+epi = ^ [ see Eqs. 03, 03]- 
Finally, we stress that the NLS Eq. 03), or the HNLS 
Eq. ([20]) . can be used in the LH (RH) regime, taking 
k , e£, and fit negative (positive) as per the discussion 
above. 

Let us now analyze Eqs. 03) an d (|20l) in more detail. 
First, measuring length, time, and the field intensity \(j>\ 2 
in units of the dispersion length = fo/|feoli initial 
pulse width to, an d Ljj/|7|, respectively, we reduce the 
NLS Eq. (f!9)) to the following dimensionless form: 



id Z (f>- -<9£</> + er|(/>| 2 = 0, 



(21) 



where s = sign(A; ') and a = sign(7). The NLS Eq. l|2T]) 
admits bright (dark) soliton solutions for scr = —1 (sa = 
+1). As is shown in Fig. [21 for our choice of parameters, 
s = +1 (i.e., k'o > 0) for 2vr x 1.76 < w < 2vr x 1.87 
GHz, while s = — 1 (i.e., fc ' < 0) for 2-zr x 1.45 < w < 
2tt x 1.76 GHz in the LH regime. Moreover, since f3 > 0, 
we have a = +1 either for a focusing dielectric, a > 0, or 
for a defocusing dielectric, a < 0, with \a//3\ < Zq/Z* 
(Zq = yVo Ao is the vacuum wave-impedance). Hence, 
for a = +1, bright (dark) solitons occur in the anomalous 
(normal) dispersion regimes, i.e., for fc ' < (fc ' > 0), 
respectively. On the other hand, a = — 1 for a < 0, 
with |a//3| > Zq/Z^ and, bright (dark) solitons occur in 
the normal (anomalous) dispersion regimes. The above 
results are summarized in Table I. Note that the presence 
of dispersion and nonlinearity in the magnetic response 
of the metamaterial allows for bright (dark) solitons in 
the anomalous (normal) dispersion regimes for defocusing 
dielectrics (see third line of Table I) . 

Next, we consider the HNLS Eq. ([20| which, by using 
the same dimensionless units as before, is expressed as, 

id z ^>- ^a^$ + cr|$| 2 $ = i8 1 d^-i<78 2 dT{\<$>\ 2( $>), (22) 

where 8\ = £fe "/(6to|fc '|), a nd <5 2 = e/(uiota). Equation 
(|22p can be used to predict ultra-short solitons in nonlin- 
ear LH metamaterials as follows. Following Ref. [2fJ, we 
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seek travelling-wave solutions of Eq. (|22| of the form, 



$(Z, T) = U{rj) exp[i(KZ - flT)} 



(23) 



where U (j]) is the unknown envelope function (assumed 
to be real), 77 = T—AZ, and the real parameters A, K and 
n denote, respectively, the inverse velocity, wavenumber 
and frequency of the travelling wave. Substituting Eq. 
([23| into Eq. (|22|) . the real and imaginary parts of the 
resulting equation respectively read: 



U 



(24) 



Silf + (A - sft - 3<5i0 2 )?7 - 3aS 2 U 2 U = 0, (25) 

where overdots denote differentiations with respect to rj. 
Notice that in the case of 5\ = 62 = Q, Eq. ([25]) is au- 
tomatically satisfied if A = s£l and the profile of "long" 
soliton pulses [governed by Eq. (|2"Tj) ] is determined by Eq. 

On the other hand, for ultra-short solitons (corre- 
sponding to 61 0, S 2 7^ 0), the system of Eqs. (|2~4")1 and 
is consistent if the following conditions hold: 



K-^n 2 - S^ 3 _A- S n- 3<5i0 2 
§ + 36 & ~~ 5[ 
aS 2 a(l + S 2 Q) _ 
~ ~6\ ~ f + 3<Sifi = 



(26) 
(27) 



where k and ^ are nonzero constants. In such a case, Eqs. 

and (|23|) are equivalent to the following equation of 
motion of the unforced and undamped Duffing oscillator, 



U + kU + isU 3 = 0. 



(28) 



For kv < 0, Eq. (|28|) possesses two exponentially lo- 
calized solutions (as special cases of its general elliptic 
function solutions), corresponding to the separatrices in 
the (U, U) phase-plane. These solutions have the form of 



a hyperbolic secant (tangent) for k < and v > (k > 
and v < 0), thus corresponding to the bright, t/ BS (dark, 
Ubs) solitons of Eq. l[2"2"jl : 

U BB (V) = (2k|/^) 1/2 sech(/Rr;), (29) 
Unsirj) = (2 K /|^|) 1/2 tanh(V^72^)- (30) 



These are ultra-short solitons of the HNLS Eq. (|22|) , valid 
even for e = O(l): since both coefficients 8\, S 2 of Eq. 
([22]) scale as e(woio) -1 : it is clear that for wo^o = C(l), or 
for soliton widths t ~ cj^ 1 , the higher-order terms can 
safely be neglected and soliton propagation is governed 
by Eq. |[2"T]). On the other hand, if cj <o = C(e), the 
higher-order terms become important and solitons gov- 
erned by the HNLS Eq. (|22|) are ultra-short, of a width 
to ~ £Uq . We stress that these solitons are approxi- 
mate solutions of Maxwell's equations, satisfying Fara- 
day's and Ampere's Laws in Eqs. (|5|) up to order 0(e 3 ). 

Finally, as concerns the condition for bright or dark 
soliton formation, namely kv < 0, we note that k de- 
pends on the free parameters K and fl (and, thus, can 
be tuned on demand), while the parameter v has the 
opposite sign from a (since S 2 > 0, while signal) = 
sign(fc ") = +1 - see Fig. [2]). This means that bright soli- 
tons are formed for k < and a = — 1 (i.e., a < with 
\a//3\ > Zq/Z^), while dark ones are formed for k > 
and a = +1 (i.e., a > 0, or a < with \a/[3\ < Z^/Z^). 

In conclusion, we used the reductive perturbation 
method to derive from Maxwell's equations a HNLS 
equation describing pulse propagation in nonlinear meta- 
materials. We studied the pertinent dispersive and non- 
linear effects, found necessary conditions for the forma- 
tion of bright or dark ultra-short solitons, as well as ap- 
proximate analytical expressions for these solutions. Fu- 
ture research may include a systematic study of the sta- 
bility and dynamics of the ultra-short solitons, both in 
the framework of the HNLS equation and, perhaps more 
importantly, in the context of Maxwell's equations. 
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